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SUMMARY 

A  systematic  study  of  threat  effectiveness  for  antiaircraft  artillery  (AAA) 
systems  requires  the  development  of  a  mathematical  model  for  the  gunner's 
tracking  response.  The  gunner  (or  tracker)  model  is  then  incorporated  into 
computer  simulation  programs  for  predicting  aircraft  attrition  with  respect 
to  specific  antiaircraft  weapon  systems.  Two  of  the  fundamental  design 
requirements  of  a  gunner  model  are  simplicity  in  model  structure  and 
accuracy  in  the  tracking  error  predictions.  A  simple  gunner  model  structure 
will  shorten  computer  simulation  execution  time.  Obviously,  accurate 
predictions  of  tracking  error  implies  model  fidelity  with  respect  to 
describing  the  gunner's  tracking  performance.  The  Luenberger  full-order 
observer  theory  has  been  applied  to  design  a  human  operator  model  (observer 
model)  for  AAA  tracker  response  which  has  been  documented  in  a  previous 
report.  This  technical  report  will  describe  the  development  of  an  anti¬ 
aircraft  gunner  model  based  on  the  Luenberger  reduced-order  observer  theory. 

It  satisfies  both  the  design  requirements  mentioned  above.  It  is  composed 
of  three  main  parts  -  a  reduced-order  observer,  a  feedback  controller,  and 
a  remnant  element. 

A  parameter  identification  program  based  on  the  least  squares  curve-fitting 
method  and  the  Gauss-Newton  gradient  algorithm  is  developed  to  systematically 
determine  the  model  parameters.  This  program  iteratively  adjusts  the 
parameter  values  to  minimize  the  error  between  the  model  prediction  of 
tracking  error  and  actual  human  tracking  data  obtained  from  manned  AAA 
simulation  experiments  conducted  at  the  Aerospace  Medical  Research  Laboratory, 
Wright-Patterson  Air  Force  Base,  Ohio.  Computer  simulation  results -of  the 
AAA  tracking  task  using  this  model  are  in  excellent  agreement  with  the 
empirical  data  for  several  aircraft  flyby  and  maneuvering  trajectories.  A 
comparison  between  this  model  and  the  optimal  control  model  by  Kleinman, 

Baron,  and  Levison  is  also  given.  This  model  is  shown  to  be  as  accurate  as 
the  optimal  control  model  in  predicting  tracking  errors.  In  addition,  the 
computer  execution  time  of  the  AAA  closed  loop  system  simulation  utilizing 
this  model  is  less  than  15  percent  of  that  using  the  optimal  control  model. 
Therefore,  this  gunner  model  car*  be  used  accurately  and  efficiently  in  the 
study  of  the  AAA  effectiveness  and  aircraft  survivability. 
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Section  I 
INTRODUCTION 

A  computer  simulation  study  of  threat  effectiveness  for  antiaircraft 
artillery  (AAA)  systems  requires  the  development  of  a  mathmetical  model  for 
the  gunner  (or  tracker  or  human  operator)  response.  The  gunner  model 
represents  the  human  operator's  control  characteristics  in  a  compensatory 
tracking  task.  It  is  then  incorporated  into  computer  engagement  simulation 
programs  [1]  for  predicting  aircraft  attrition  with  respect  to  specific 
antiaircraft  weapon  systems.  Two  of  the  fundamental  design  requirements  of 
a  gunner  model  are  simplicity  in  model  structure  and  accuracy  in  the 
tracking  error  predictions.  A  simple  gunner  model  structure  will  shorten 
computer  simulation  execution  time.  Obviously,  accurate  predictions  of 
tracking  error  implies  model  fidelity  with  respect  to  describing  the 
gunner's  tracking  performance.  Then,  the  manned  threat  quantification  in 
the  threat  analysis  will  be  reliable. 

An  antiaircraft  gunner  model  based  on  the  Luenberger  full-order  observer 
theory  [2]  was  developed  and  documented  in  a  previous  report  [3],  In  this 
report,  the  Luenberger  reduced-order  observer  theory  [4],  [5]  is  applied  to 
develop  a  tracker  model  for  AAA  compensatory  tracking  task.  It  satisfies 
both  the  design  requirements  mentioned  above.  The  structure  of  the  model 
is  simple  and  its  predictions  of  tracking  errors  are  accurate.  It  is 
composed  of  three  main  parts  -  a  reduced-order  observer,  a  feedback  con¬ 
troller,  and  a  remnant  element.  An  observer  is  itself  a  dynamic  system 
whose  output  can  be  used  as  an  estimate  of  the  state  of  a  given  system.  A 
reduced-order  observer  has  dynamic  order  less  than  the  observed  system  and 
provides  an  estimate  of  those  state  components  which  are  not  available  for 
direct  measurement.  The  structure  of  a  reduced-order  observer  is  simple 
and  its  design  is  easy.  The  idea-  of  using  a  reduced-order  observer  in  the 
tracker  model  design  is  to  obtain  an  appropriate  estimate  of  the  state 
components  (which  are  not  directly  measurable)  of  the  gunsight  system  and 
the  target  motion.  When  the  gunner  (a  human  operator)  observes  the  tracking 
error  from  the  visual  display,  he  not  only  obtains  the  tracking  error 
information,  but  also  has  a  certain  understanding  or  knowledge  about  other 
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variables  (state  components)  of  the  overall  system.  It  is  one  of  the  main 
differences  between  human  tracking  and  machine  tracking.  A  human  operator 
(gunner)  can  always  realize  more  information  about  the  system  than  what  is 
on  the  display.  This  fact  is  represented  by  a  reduced-order  observer  in 
the  gunner  model.  The  estimated  state  components  from  the  reduced-order 
observer  and  the  observed  state  components  from  the  display  are  then  used 
to  implement  a  linear  feedback  controller  which  represents  the  gunner's 
control  function  in  the  compensatory  tracking  task.  The  effects  of  all  the 
randomness  sources  due  to  human  psychophysical  limitations  and  of  modelling 
errors  are  lumped  into  one  random  remnant  element  in  this  model  design. 
Another  important  feature  of  this  model  is  that  its  parameters  can  be 
determined  systematically  instead  of  by  trial-and-error.  A  parameter 
identification  program  based  on  the  least  squares  curve-fitting  method  [6] 
and  the  Gauss-Newton  gradient  algorithm  [7]  is  developed  for  this  purpose. 
This  prof., -am  iteratively  adjusts  the  parameter  values  to  minimize  the  least 
squares  error  between  the  model  prediction  of  tracking  error  and  actual 
human  tracking  data  obtained  from  manned  AAA  simulation  experiments  conducted 
at  the  Aerospace  Medical  Research  Laboratory,  WPAFB,  Ohio.  Thus  it  provides 
a  convenient  procedure  for  model  validation.  In  addition,  a  computer 
simulation  program  ROOMS  (Reduced-Order  Observer  Model  Simulation)  is 
developed  with  the  designed  model  describing  the  gunner's  response  for  a 
given  AAA  tracking  task.  The  program  provides  time  functions  of  the 
ensemble  mean  and  standard  deviation  for  the  model's  tracking  error  pre¬ 
dictions  (azimuth  and  elevation) .  Computer  simulation  results  are  in 
excellent  agreement  with  the  empirical  data.  Furthermore,  this  model  is  a 
predictive  model  in  the  sense  that  it  can  be  used  to  predict  tracking 
errors  of  an  AAA  system  for  various  flyby  and  maneuvering  trajectories  with 
similar  frequency  bandwidths. 

A  comparison  between  this  model  and  the  optimal  control  model  [8],  [9],  and 
[10]  (by  Kleinman,  Baron,  Levison)  is  also  given.  It  can  be  shown  that  the 
model  based  on  observer  theory  is  as  accurate  as  the  optimal  control  model 
in  predicting  tracking  errors.  In  addition,  the  computer  execution  time  of 
the  AAA  closed  loop  system  simulation  utilizing  this  model  is  less  than 
15  percent  of  that  using  the  optimal  control  model.  This  is  a  primary 
advantage  of  a  model  with  simple  structure. 
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The  design  of  the  antiaircraft  gunner  model  based  on  the  reduced-order 
observer  theory  is  described  in  detail  in  Section  II.  Section  III  gives 
the  model  validation  method  and  computer  simulation  results.  The  conclusion 
is  given  in  Section  IV. 
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Section  II 

REDUCED-ORDER  OBSERVER  MODEL 


In  [3],  the  azimuth  or  elevation  gunsight  dynamics,  rate-aided  control 
dynamics  and  the  target  motion  of  an  antiaircraft  artillery  (AAA)  gun 
system  have  been  represented  by  the  following  state  space  equation. 

x  -  Ax  +  Bu  +  F0T  (1) 
where  x  is  the  state  vector  with  two  components, 


x(t) 


Xl(t) 

X2(t) 


eT(t) 

e.f(t) 


eT(t)  is  the  tracking  error,  i.e.,  the  difference  between  the  target 
position  angle  8  and  the  gunsight  line  angle  0  .  0  is  the  target  angle 

^  **  o  ^ 

rate.  The  scalars  u  and  0T  in  Eq  (1)  denote  the  tracker’s  control  output 
and  target  angle  acceleration  respectively.  The  matrices  A,  B,  and  F  are 


The  tracking  error  eT  on  the  visual  display  is  observed  by  the  gunner  and 
is  expressed  in  the  measurement  equation: 

y  =  Cx  (2) 
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where  y  is  the  observed  tracking  error  and  C  is  a  row  vector  [1  0].  Next, 
Equations  (1)  and  (2)  will  be  used  to  develop  the  reduced-order  observer 
model. 

The  structure  of  the  reduced-order  observer  model  is  shown  in  Figure  1.  It 
consists  of  three  main  elements:  a  reduced-order  observer,  a  controller, 
and  a  remnant  element.  The  reduced-order  observer  processes  the  tracker’s 
observation  from  the  visual  display  to  provide  an  estimate  of  those  state 
components  of  the  AAA  system  which  are  not  directly  measurable.  It  will  be 
shown  that  the  system  equation  (1)  is  a  second  order  system,  but  the  reduced- 
order  observer  is  only  a  first  order  system,  since  some  components  of  the 
state  vector  as  given  by  the  system  outputs  are  already  available  by  direct 
measurement.  The  estimation  of  those  measurable  state  components  is  nof. 
necessary  and  will  cause  a  certain  degree  of  redundancy.  The  use  of  a 
reduced-order  observer  eliminates  this  redundancy  and  provides  an  approximate 
estimation  of  the  state  components  which  can’t  be  measured  directly.  The 
controller  represents  the  gunner’s  tracking  function  by  a  state  variable 
linear  feedback  control  law.  The  observer  c.d  the  controller  consists  of 
the  deterministic  part  cf  the  gunner  model.  The  effects  of  the  various 
randomness  sources  in  the  AAA  man-machine  closed  loop  system  and  of  the 
modelling  errors  are  lumped  into  one  element  called  remnant,  which  is  the 
stochastic  part  of  the  gunner  model.  These  randomness  sources  include  the 
modelling  error,  the  observation  error,  the  neuromotor  noise,  etc.  Mathema¬ 
tical  equations  of  this  model  are  given  below. 

A.  Redueed-Oraer  Observer  Design 

System  equations  (1)  and  (2)  are  used  in  the  design  of  the  first  element 
(reduced-order  observer)  of  the  gunner  model.  However,  the  gunner  does 
not  have  precise  information  about  the  target  dynamics,  so  the  term  repre¬ 
senting  target  acceleration,  0T>  in  Equation  (1)  will  not  be  included  in 
the  design  of  the  observer  equation.  The  effect  on  the  tracking  error  due 
to  eliminating  the  0^  term  will  be  included  in  the  remnant  element.  Now 
from  Equation  (2),  y  =  Cx  =  xi ,  the  cracking  error  is  available  from  direct 
observation.  Thus,  it  >.•>'  only  necessary  to  estimate  the  second  component 
x 2  of  the  state  vector  x  in  order  to  implement  a  state  variable  feedback 
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Figure  1.  Block  Diagram  of  the  Structure  of  the  Reduced-Order  Observer  Model 


control  law.  In  the  following,  the  Luenberger  reduced-order  observer 
theory  [4]  and  [5]  is  used  to  design  the  gunner  model.  First,  with  the  6^ 
term  eliminated,  Equation  (1)  can  be  rewritten  as 

X1  =  allxl  +  ai  2X2  +  biu  (3) 

x2  =  a21xl  +  a22x2  +  b2u  (4) 

Since  the  first  state  component  xi  is  measurable,  i.e.,  y  =  xj ,  Equation  (3) 
can  also  be  expressed  by 

y  *  auy  +  a12x2  +  b2u 


or  equivalently. 


y  -  any  =  ai2x2  +  biu 


(5) 


Let  us  introduce  a  new  variable  y'  *  y  -  aj iy,  then  Equations  (4)  and  (5) 
can  be  expressed  by 


x2  =  a22x2  a2iy  b2U  (6) 

and 

y'  ■  ai  2x2  +  bi,u  (7) 

Now  Equation  (6)  is  the  reduced-order  system  dynamic  equation  with  measure¬ 
ment  data  obtained  by  Equation  (7).  Note  that  Equations  (6)  and  (7)  are  a 
first  order  system  with  one  measurement  equation.  An  observer  which  gives 
an  estimate  £2  of  x2  can  be  easily  designed  as  shown  by  the  following 
equation. 


9  =  a99Xo 


*  (a22  -  ka12)  x2  +  ky  +  (a2I  -  kan)  y  +  (b2  -  M^)  u£ 


(8) 
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where  and  are  the  elements  of  matrices  A  and  B  in  Equation  (1) ,  the 
scalar  k  is  the  observer  gain,  y  and  y  are  the  observed  tracking  error  and 
error  rate  respectively,  and  u^  is  the  linear  feedback  control  law  (the 
controller)  with  the  form: 

“c  =  -  **][*,] 

where  the  feedback  control  gains  yj  and  Y2  are  two  constants  to  be  determined 
later.  Note  that  the  state  feedback  is  composed  of  y  (the  observed  variable 
which  is  xj)  and  $2  (the  estimated  state  of  X2).  It  can  be  shown  that  the 
system  (1)  and  (2)  is  completely  observable.  (The  definition  of  observability 
and  the  conditions  of  a  system  to  be  observable  can  be  found  in  [11],  Then, 
by  the  observer  theory,  there  always  exists  an  observer  gain  k  to  make  the 
eigenvalue  of  the  observer  (Equation  (8))  negative.  Thus,  the  output  of  the 
observer  will  be  a  good  estimation  to  the  state  of  the  observed  system. 

This  shows  the  existence  of  proper  observer  gain  k  in  Equation  (8) .  Actually, 
the  value  of  observer  gain  k  is  determined  by  a  curve-fitting  identification 
program.  The  required  differentiation  of  y  in  Equation  (8)  can  be  avoided 
by  introducing  the  following  variable: 

z(t)  =  k2  -  ky(t)  (9) 

Hence  the  observer  dynamics  can  be  represented  by 

z  =  (a22  -  kai2)  z  +  (a22  -  kai2)  ky  +  (a2i  -  kan)  y  + 

(b2  -  kbj)  uc  (10) 

Next,  the  actual  output  of  this  model  is  expressed  as  the  sum  of  the  output 

u  of  the  controller  and  the  remnant  element  v. 
c 

u  =  u  +  v 
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where  the  remnant  term  v(t)  is  modeled  as  a  white  noise  and  its  statistical 
properties  are  selected  to  be 

E  [v(t)]  =  0  for  all  t 

E  fv(t)  v(t)]  =  q(t)  6 (t  -  t)  for  all  t  and  x  (12) 

where  E  is  the  expectation  operator  and  the  covariance  function  q(t)  is  a 
scalar  function  of  time  and  will  be  described  later.  6  is  the  Dirac  delta 
function. 


The  gunner  model  equations  of  the  observer,  the  controller,  and  the  remnant 
have  been  derived.  These  equations  are  combined  with  system  equations  (1) 
and  (2)  to  obtain  the  mathematical  model  of  the  closed  loop  AAA  system. 
Since  x^  =  y,  Equations  (1)  and  (10)  can  be  rewritten  as  follows 


y  =  any  +  ai2x2  +  biU  +  f  10  ' 

T  (13) 

•  ** 
x2  =  a21y  +  a22x2  +  b2u  +  f20T 

z  =  (a22  -  ka12)  z  +  (a22  -  ka12)  ky  +  (a21  -  kan)  y  + 

(b2  -  kbi)  u£ 

u  =  u  +  v 
c 

uc  -  -  [«  «][*J 
By  introducing  new  variables: 
x3  =  x2  -  ky 
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and 


e  =  x3  -  z 


(14) 


Equation  (13)  can  be  rewritten  as 


y  =  [an  +  a12k  r  bi (yx  +  ky2)J  y  +  (a12  -  bjY2)x3  + 


biY2e  +  fi0T  +  biv 


x3  =  [(a22  “  kan)  k  +  (a2i  “  kan)  -  (b2  -  kbi)  (yi  +  ky2)J  y  + 

[a22  -  kai2  -  (b2  -  kbj)  Y2J  x3  +  (b2  “  kbi)  Y2e  + 

(b2  -  kbx)  v  +  (f2  -  kfi)  0T 

•  • 

e  =  (a22  -  ka12)  e  +  (f2  -  kfi)  0T  +  (b2  -  kb:)  v 


Or  equivalently, 


X  =  AxX  +  Fx0  +  Dxv 


(15) 


where  X  is  the  state  vector  of  the  overall  system  with  components: 


y 

y 

X  = 

x3 

- 

x2  "  ky 

e 

x3  -  z 

An  Fn  and  Dj  are  matrices  defined  as  follows: 
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Ai  = 


all  +  a12k  "  bl  (Yl  +  k  Y2) 


(a22  -  ka12)  k  +  a21  -  kan 
(b2  -  kbj) (yi  +  ky2) 


Fl 


fl 

f2  -  kfj 

«2  -  kfl 


a12  “  bl  Y2 


a22  “  ka12  - 

(b2  -  kbx)Y2 


Dj  = 


bl 

b2  -  kbj 
b2  -  kbi 


bi  Y2 


(b2  -  kbx)y2 


a22  -  ka12 


In  order  to  determine  the  parameters  k,  yj,  y2>  and  the  covariance  function 
q(t)  in  Equation  (12),  we  use  the  following  equations.  Letting  the 
expectation  value  of  X  be  X  then  we  have, 


X  -  AiX  +  FX0T 


(16) 


and  the  covariance  matrix  of  X(t)  is  P(t)  =  E  [( x( t )  -  x(t))  (x(t)  -  x( t) ) T] ; 
then  it  can  be  shown  in  [  11  ]  that  the  covariance  matrix  is  governed  by 


P  =  AjP  +  PAj1  +  Dj  q(t)  Di1 


(17) 


Equations  (16)  and  (17)  will  be  used  in  Section  III  to  determine  the  values 
of  parameters  in  the  gunner  model. 
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B.  Frequency  Domain  Analysis 

A  detailed  study  of  the  frequency  domain  responses  of  some  key  time  domain 

variables  of  the  AAA  tracking  system  has  been  done.  This  will  provide  us 

information  to  design  the  covariance  function  q(t)  of  the  remnant  element 

(see  Equation  (12)).  The  power  spectral  density  functions  (PSD)  of  the 

•  ** 

target  angle  8^,  angle  rate  8^  and  angle  acceleration  8^  were  generated  for 
the  four  flyby  and  maneuvering  trajectories  [3],  These  PSD  functions  were 
generated  for  both  azimuth  and  elevation  components  of  target  motion 
variables.  It  was  found  that  the  significant  parts  of  all  these  PSD 
functions  are  less  than  0,5  Hz.  In  addition,  similar  PSD  functions  were 
generated  for  the  sample  ensemble  mean  e^  of  tracking  errors  which  were 
obtained  by  averaging  empirical  data  of  sixteen  runs.  The  empirical  data 
was  generated  from  manned  AAA  simulation  experiments  conducted  at  the 
Aerospace  Medical  Research  Laboratory,  WPAFB.  There  is  a  significant  con¬ 
sistency  between  the  PSD's  of  target  angle  acceleration  0T  and  the  PSD's  of 
the  corresponding  sample  ensemble  mean  e^,  of  the  tracking  error.  These 
results  are  shown  in  Figures  2  through  9.  There  are  two  curves  on  each  of 
these  figures.  The  solid  curve  denotes  the  PSD  function  of  the  empirical 
mean  of  tracking  error.  The  dashed  curve  describes  the  corresponding  PSD 
function  of  the  target  angle  acceleration.  Figures  2  through  5  show  the 
comparison  of  PSD  functions  of  the  azimuth  case  for  four  target  trajectories, 
respectively.  Figures  6  through  9  show  the  similar  results  for  the  elevation 
case.  The  low  frequency  parts  of  the  PSD  functions  of  target  accelerations 
match  well  the  low  frequency  parts  of  the  PSD  functions  of  empirical  means  • 
of  corresponding  tracking  errors.  The  high  frequency  parts  of  the  PSD 
functions  of  empirical  mean  tracking  errors  in  these  figures  are  the  results 
due  to  small  sample  size  of  experimental  runs.  Figure  8  shows  that  these 
two  PSD  functions  do  not  match.  After  further  investigation,  it  was  found 
that  the  data  tape  containing  the  information  of  target  elevation  acceler¬ 
ation  was  not  properly  generated.  In  this  frequency  domain  analysis,  it  is 
concluded  that  the  tracking  error  in  AAA  tracking  task  depends  on  the 
target  acceleration.  In  other  words,  when  the  target  aircraft  makes  a 
maneuvering  flight  (i.e.,  the  corresponding  target  acceleration  increases), 
the  tracking  difficulty  increases  and  hence  the  tracking  error  increases. 

This  observation  will  be  used  in  the  next  section  to  design  the  covariance 
function  of  remnant  element. 
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EMPIRICAL  mean  of 
tracking  error 


Trajectory  2  (Azimuth) 


Trajectory  3  (Azimuth) 


empirical  mean  of 
TRACKING  ERROR 


Comparison  of  PSD  Functions  Trajectory  1  (Elevation) 


Figure  7.  Comparison  of  PSD  Functions  Trajectory  2  (Elevation) 


EMPIRICAL  MEAN  OF 


Figure  8.  Comparison  of  PSD  Functions  Trajectory  3  (Elevation) 


Section  III 
MODEL  VALIDATION 

A.  Time  Domain  Curve-Fitting  Identification  Method 

The  design  of  a  gunner  model  included  two  main  procedures;  the  design  of 
the  structure  of  the  model  and  the  determination  of  the  values  of  the  model 
parameters.  The  structure  of  the  gunner  model  based  on  observer  theory  has 
been  described  in  Section  II.  In  this  section,  the  parameters  of  the  model 
will  be  determined  through  an  identification  program.  Figure  10  shows  the 
parameter  determination  procedure.  The  manned  AAA  simulator  built  at  the 
Aerospace  Medical  Research  Laboratory,  Wright-Patterson  AFB,  Ohio,  provides 
the  empirical  tracking  error  data,  e^t),  along  with  the  simulated  target 
trajectory  8^  as  the  input.  e^,(t)  represents  the  sample  ensemble  mean  of 
the  tracking  error  of  sixteen  simulation  runs.  Now  in  Equation  (16),  the 
first  component  of  X  is  the  expectation  value  of  the  tracking  error,  i.e., 

_  I 

the  model  prediction  e^  of  the  ensemble  mean  of  the  tracking  error  in  the 
AAA  closed  loop  system.  Ae  denotes  the  difference  between  the  empirical 

_  _  i 

data  eT  and  the  model  prediction  eT«  The  criterion  function  evaluates 
the  "goodness  of  fit"  between  the  model  prediction  and  the  empirical  data. 

In  this  study,  the  criterion  function  is  selected  to  be 

J  =/0f  (eT(t)  -  eT  (t,  y  »  Y2>  k))2dt 

_  f 

where  t^  is  the  tracking  duration  (equal  to  45  seconds) ,  e^  is  the  model 

prediction  of  the  tracking  error  which  is  a  function  of  time  and  the  unknown 

_» 

parameters.  The  explicit  form  of  the  function  eT  in  terms  of  time  t  and 

model  parameters  y  ,  y  and  k  has  been  derived  and  is  included  in  Appendix  A. 
1  2 

The  criterion  function  J  will  be  minimized  with  respect  to  these  unknown 
parameters.  The  Gauss  Newton  gradient  algorithm  was  used  as  the  parameter 
adjustment  algorithm  to  iteratively  determine  the  values  of  model  parameters. 
The  values  of  parameters  obtained  through  this  identification  program  are: 
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Figure  10.  Block  Diagram  of  the  Parameter  Identification  Procedure 


Azimuth  Tracking 

Elevation  Tracking 

Observer  gain  k 

2.94 

3.02 

Controller  gains  y^ 

-2.87 

-3.01 

y 

-1.00 

-1.00 

2 

Once  the  observer  gain  and  the  controller  gains  are  determined,  the  system 
matrix  Ai  is  known.  Next  the  model  prediction  of  the  standard  deviation  of 
the  tracking  error  is  considered.  It  can  be  shown  that  the  square  root  of 
the  first  diagonal  element  Pn  of  the  covariance  matrix  P(t)  in  Equation  (17) 
is  the  standard  deviation  of  the  tracking  error.  The  solution  of  Equation 
(17)  is 


P(t)  =  *(t,  tQ)  P(tQ)  *T(t,  to)  + 

f\  <K t,  T)Diq(x)D1V(t,  T)dT  (18) 

where  4>(t,  tQ)  is  the  state  transition  matrix  defined  as 
<t>(t,  tQ)  *  eAl^t“to-) 

and  q(t)  is  the  covariance  matrix  of  the  random  remnant.  It  has  been  found 

from  the  frequency  domain  analysis  results  that  the  tracking  error  depends 

•  • 

on  the  target  trajectory  dynamics,  especially  the  target  acceleration  0^. 
Therefore,  it  is  assumed  that  the  covariance  function  q(t)  of  the  remnant 
is  a  function  of  the  target  dynamics  as  follows: 

q(t)  =  dj  +  ct20T2(t)  +  a3eT2(t)  (19) 

where  a  ,  a  ,  and  a  are  three  nonnegative  constants  to  be  determined  later 
1  2  3 
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and  0_  and  0_  are  estimated  target  angle  rate  and  acceleration,  respectively. 
T  1  j  .'i 

The  reason  that  only  estimated  values  0^  and  0^  were  used  in  Equation  (19) 
is  that  the  gunner  does  not  have  precise  information  about  0^  and  0  (i.e., 

target  uncertainty).  These  estimated  quantities  can  be  obtained  as  follows: 
Equation  (9)  can  be  rewritten  by 


0T  =  z(t)  +  ky(t) 


where  z(t)  is  the  output  of  the  reduced-order  observer,  k  is  the  observer 

gain,  and  y(t)  is  the  observed  tracking  error.  Since  all  these  quantities 

*  g 

are  known,  0  can  be  computed.  Next,  by  the  first-order  approximation,  0^ 

can  be  found  from, 


9  T^tk) 


A  A 

0  T  (tk)-  6T(tk-1) 
At 


where  tk  *  k&t,  and  At  *  0.04  sec  is  the  sampling  period.  Then  by  minimizing 

the  following  cost  function  J*,  the  parameters  a  ,  a  ,  and  a  can  be 

1  2  3 

determined, 

J’  *  f  J  [S  (t)-  S  (t,  a  ,  a  ,  a  )]2dt 
J  0  1  2  12  3 


where  is  the  sample  ensemble  standard  deviation  of  empirical  tracking 

error  of  sixteen  runs  and  S  is  the  model  prediction  of  the  ensemble  standard 

2 

deviation  (an  explicit  function  of  time  and  the  parameters  a^,  a^,  and  cO . 

Note  that  the  curve-fitting  of  the  standard  deviation  of  tracking  error  is 

also  done  in  the  time  domain  instead  of  in  the  frequency  domain  as  in  [3]. 

The  derivation  of  the  function  S  from  Equation  (18)  is  included  in 

2 

Appendix  B.  The  results  of  the  parameter  determination  program  in  minimizing 
J'  are, 
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Azimuth 

Elevation 

a 

.0496 

.0032 

1 

a 

.0024 

.00047 

2 

a 

.103 

.259 

3 

B.  Computer  Simulation  Results 

The  numerical  values  of  the  parameters  of  this  gunner  model  were  determined 
in  the  previous  subsection  with  respect  to  the  gunsight  dynamic  system 
(Equation  (1))  and  a  deterministic  target  trajectory.  The  gunner  model  is 
now  ready  to  be  used  for  computer  simulation.  A  computer  program  called 
ROOMS  simulates  the  AAA  system  with  this  model  representing  the  gunner 
response.  The  input  to  this  program  is  the  target  motion  trajectory.  The 
outputs  are  the  model  predictions  of  the  ensemble  mean  and  standard  deviation 
of  the  tracking  error.  These  results  are  plotted  in  Figures  11  through  26. 
There  are  two  curves  in  each  of  these  Figures.  The  solid  curve  denotes  the 
empirical  tracking  data.  The  corresponding  model  prediction  is  denoted  by 
the  dotted  curve. 

Figures  11  through  14  give  comparisons  between  the  empirical  data  and  the 
model  prediction  of  the  ensemble  mean  of  azimuth  tracking  errors  for  four 
target  trajectories.  Obviously,  the  model  predictions  match  the  empirical 
data  very  well  for  both  flyby  and  maneuvering  trajectories.  It  should  be 
emphasized  that  these  results  are  obtained  by  using  the  designed  gunner 
model  with  the  same  set  of  model  parameters  (determined  in  the  previous 
subsection)  applied  to  various  trajectories.  In  other  words,  it  verifies 
that  with  the  same  set  of  parameters,  the  gunner  model  based  on  the  reduced- 
order  observer  theory  can  give  accurate  predictions  of  tracking  errors  in 
AAA  weapon  systems  for  various  target  trajectories  with  similar  frequency 
bandwidths.  Therefore,  it  is  a  predictive  model.  Next,  Figures  15  through 
18  give  comparisons  between  the  empirical  data  and  the  model  prediction  of 
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Figure  11.  Mean  Tracking  Error  Trajectory  1 
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Figure  13.  Mean  Tracking  Error  Trajectory  3 


Figure  14.  Mean  Tr. 
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Figure  15.  Standard  Deviation  of  Tracking  Error  Trajectory  1 
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Figure  16.  Standard  Deviation  of  Tracking  Error  Trajectory  2 


Figure  17.  Standard  Deviation  of  Tracking  Error  Trajectory  3 
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MODEL  PREDICTION 


Figure  21.  Mean  Tracking  Error  Trajectory  3 
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Figure  23.  Standard  Deviation  of  Tracking  Error  Trajectory  1 
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Figurt  25.  Standard  Deviation  of  Tracking  Error  Trajectory  3 


the  ensemble  standard  deviation  of  azimuth  tracking  errors  for  four 
trajectories.  Again,  excellent  matching  between  the  two  curves  are  shown 
in  each  of  these  four  figures.  Similar  results  of  the  elevation  case  are 
shown  in  Figures  19  through  26  for  the  four  trajectories.  All  these 
results  indicate  that  the  gunner  model  is  able  to  represent  the  character¬ 
istics  of  the  gunner  response  in  the  AAA  compensatory  tracking  task. 

C.  Comparison 

A  comparison  of  the  model  prediction  accuracy  between  this  gunner  model  and 
the  optimal  control  model  has  been  done  for  several  target  trajectories. 

All  the  results  show  that  both  models  give  accurate  predictions  of  tracking 
errors.  Some  typical  results  are  shown  in  Figures  27  through  30.  It  is 
obvious  that  the  gunner  model  developed  by  the  authors  can  predict  the 
tracking  errors  as  accurately  as  those  obtained  by  the  optimal  control 
model.  However,  the  computer  execution  time  for  simulating  the  AAA  gun 
system  using  the  gunner  model  is  less  than  15%  of  that  used  by  the  optimal 
control  model.  It  is  a  primary  advantage  of  a  model  with  simple  structure. 
The  following  table  shows  some  highlights  of  the  gunner  model  over  the 
optimal  control  model. 


Gunner  Model 
(Observer  Theory) 

Optimal  Control 

Model 

Model  Structure 
Complexity 

Low 

High 

Computer  Simulation 

Time  (Seconds) 

5.48 

37.02 

Model  Predictions 
of  Tracking  Errors 

Good 

Good 

Model  Validation 

Parameter 

Identification 

Program 

Trial-and-Error 

Tuning 
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Section  IV 
% 

CONCLUSION 

The  Luenberger  reduced-order  observer  theory  has  been  applied  to  design  an 
antiaircraft  gunner  model  which  is  composed  of  a  reduced-order  observer,  a 
state  variable  feedback  controller  and  a  remnant  element.  The  highlights 
of  this  model  are  simple  structure  and  accurate  model  predictions.  The  key 
design  requirement  is  to  make  the  model  structure  simple  so  that  it  will 
shorten  computer  simulation  time.  It  has  also  been  shown  in  Figures  11 
through  30  that  this  model  can  predict  the  tracking  errors  accurately.  In 
addition,  parameter  identification  program  based  on  the  least  squares 
curve-fitting  method  and  the  Gauss  Newton  algorithm  has  been  developed.  It 
provides  a  systematic  procedure  to  determine  the  numerical  values  of  the 
model  parameters.  This  gunner  model  has  been  used  to  study  the  AAA 
effectiveness  of  several  foreign  air  defense  weapon  systems  at  the  Aerospace 
Medical  Research  Laboratory,  Wright-Patterson  AFB.  All  the  results  show 
that  it  is  an  accurate  and  efficient  antiaircraft  gunner  model.  Therefore, 
this  model  can  be  used  to  study  AAA  effectiveness  and  aircraft  survivability. 
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•  APPENDIX  A 


Derivation  of  Model  Prediction  of  Mean  Tracking  Error 

The  model  prediction  1^,  of  the  ensemble  mean  of  the  tracking  error  as  an 

explicit  function  of  time  t  and  parameters  y  ,  y  ,  and  k  is  derived  in  this 

1  2 

Appendix.  First  rewrite  Eq.  (16) 

X  =  A  X  +  F  0_ 

1  1  T 

It  can  be  shown  that  the  solution  X  of  Eq.  (16)  is  given  by 


X(t)  -  f  4>(t  -  t)f  eT(T)dt 

J  -00  1  L 

where  |  is  a  3  x  3  matrix  and  equals: 

<|>  <l>  <j> 

11  12  13 

<f>  <t>  <t> 

21  22  23 

4>  4> 

31  32  33 


The  first  component  Xj (t)  represents  the  ensemble  mean  of  the  tracking 
error.  From  Eqs.  (1)  and  (15), 
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iT(t)  can  be  expressed  by 


hit  -  T)  +  <|) 

(t  -  T)1 

J  —00 

L  12 

13  J 

&T(x)dT 


Let  G(s)  be  a  3  x  3  matrix  which  represents  the  Laplace  transform  of  the 
transition  matrix  <Kt);  then  it  can  be  shown  that 


G(s)  =  (si  -  Ax) 


.j  Adj(sl  -  Aj) 
det(sl  -  Aj) 


where  s  is  the  variable  of  the  complex  plane  and  I  is  a  3  x  3  identity 
matrix,  det  and  Adj  denote  the  determinant  and  adjoint  of  a  matrix 
respectively.  Let  G^s)  be  the  ijth  element  of  the  matrix  G(s) ;  then 


G  (s)  = 

12 


1  +  Y2 
S2  -  Y,8 


Gn(s) 


s2  +  (k  -  yn)s  -  ky. 

*  X 


For  simplicity,  select  y2  “  -1*  then  G^s)  "  0  and 


G  (s)  *  — - 

13  s2  +  (k  -  YjJs  -  kYj 
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Taking  the  inverse  Laplace  transform  of  G13(s),  we  have 


*  (t) 

13 


y i -kt 

1  —  p 


e  1  -  e 

y  +  k 

1 


»  ,  C  t  eYl(t  "  T)  -  e"kCt  "  t)  •• 

eT(t)  -  J  S- - — jrf -  eTMdT 

*  —oo  1  1 


This  is  the  explicit  function  of  e,p  in  terms  of  the  unknown  parameters. 

* 

The  derivicives  of  ex(t)  with  respect  to  parameters  k  and  Yj  are  expressed 
in  the  following 


3eT(t) 


•'-00  (Yx  +  k)2 


eT(T)dT 


3ex(t)  -t  (Yj  +  k)  (t 


ri 


-  -r)eYl<t  -  T)  -  (ey1(t  -  T)  -  e-k(t  -  O) 
(Yx  +  k)2 


0T(T)dx 


These  quantities  are  used  in  the  mean  curve-fitting  parameter  identification 
program. 


ki'M  1 


.4  A 


APPENDIX  B 


Derivation  of  Model  Prediction  of  Standard  Deviation  of  Tracking  Error 


The  function  S2  which  is  the  square  of  the  standard  deviation  of  the 

tracking  error  is  derived  in  terms  of  a  ,  a  ,  and  a  .  Equation  (18)  is 

12  3 

rewritten  here: 


P(t)  »  <j>(t,to)P(to)<}>  (t,tQ)  + 


f  f(t,t)D  q(t)D  V(t,T)dr 

J  t  11 


where 


q(t)  *  a  +  «  0T2(O  +  «  L2(t), 


'  bl  ■ 

'-I  * 

D1  ' 

b2  -  kb, 

B 

k 

b2  -  kb 

k 

• 

and  let 


Q(t)  =  D  q(r)D  T 


1 

-k 

-k 

-k 

k2 

k2 

-k 

k2 

k2 

q(x) 


Since  the  first  element  p  in  the  diagonal  of  matrix  P(t)  is  the  square  of 
the  standard  deviation  of  the  tracking  error  and  denoted  by  S2,  then. 


s?  =  *„<*> =  A  *n2  <t>to)Pii(to)  +  si  A  A 


i=l 


-o  j-1  k=l 


4>lk(t_T>qkj 

where  q,  .  denotes  the  kjth  element  of  the  matrix  Q.  As  shown  in  Appendix  A, 

the  parameter  y  is  selected  to  equal  -1.  Therefore  $  (t)  =  0.  Then 
2  12 

S  (t)  =  <j>  2(t,t  )p  (t  )  +  <j>  2(t,t  )p  (t  )  + 

2  11  0  11  0  13  0  33  o 

/[ 4>  2(t--r)  +  <j>  2(t-x)k2  -  2  <J>  (t-x)  <)>  (t-x)k]  q  (x)dx 

t  11  13  11  13 


where  <f>  (t)  has  been  computed  in  Appendix  A  and  (f>  (t)  is  obtained  in  the 

13  11 

following.  From  Appendix  A,  it  can  be  shown  that 


G  (s)  =  — - - 

11  s  -  y 


So  we  have 


Y,t 

4*i  1  (t;  =  e  a 


Then 


S  (t)  =  e2Yl(t-to)p  (t  )  +  - 


Y! (t-t0)  -k(t-t0) 


-  e 


11  0 


Yj  +  k 


P  (t  )  + 

33  0 
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Y 


•  ^  +  ct^  §T2  (t)  +  0’T2  (x)j  dr 


This  is  the  explicit  function  of  S  in  terms  of  a  ,  a  ,  and  a  . 

2  12  3 


£,(t)  ‘  e  ‘  °  pn  C°  +  l” - - 


£  (t)  * 

2 


Y  t  ^ -kt 

,1  -  k  e  ■  -  -e _ 

Yx  +  k 


(t)  *  f  (t)  +  f  f  (t-x)  (a  +  a  0  2(x)  +  a 

!  1  *  2  \  1  2  T 


3'0t  (X)  ]  dx 


Now  the  derivations  of  S  (t)  with  respect  to  a  ,  a  ,  and  a  are 

2  12  3 


as2(t)  r t 

—  -  J  t  vf-t)dT 

1  0 


as2(t) 


/t  i 

t  f2(t-x)0T2(x)dx 


3S2(t) 


_  f  (t-x)e  2(x)dx 


These  quantities  are  used  in  the  standard  deviation  curve-fitting  parameter 
identification  program. 


APPENDIX^ C 


Flow  Chart  and  Program  Listing  of  Time  Domain  Curve-Fitting  Program 

A.  Flow  Chart  of  Parameter  Identification  for  Parameters  Associated  with 
Mean  Tracking  Error  Equation. 
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Comments:  Execution  of  Program  Fit. 


1.  Compile  Program 

2.  Input:  Tape  1  contains 

a.  Azimuth  and  elevation  trajectory  information 

b.  Empirical  data  -  mean  tracking  error 

3.  Load  Library  Routines: 

a.  VCONVO  -  convolution  integral  computation  from  International 
Mathematical  &  Statistical  Library  (IMSL) . 

b.  GMINV  -  matrix  inverse  routine  from  D.  L.  Kleinman  Library 

[12]. 

4.  Execute 

5.  Output:  Parameter  values  [k,  Yj], 
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1 


PROGRAM  FIT (INPUT , OUTPUT, TA^Ei) 

C0MM0N/0AR/LIM1,LIN2,£E,IP,H,N,0EL,N1,  N2  ’ 

C  OH  HON  /M  A  T /  F2  ( i‘J  2'M )  ,  A  ( 2)  ,  AA  (  2 )  ,  B  ( 10  2<  ,  2)  ,  0  { 2)  ,  F3  ( A  096) 

DATA  N,H,Ia»Nl»  (2/1024,  .Ji*, 2, 100, 1024/  "  < 

C  CURVE  FIT  ROUTINE  TO  DETERMINE  IP  =2  PARAMETERS  KGAIN  S,  GAMMA1  i 

ffo"  id  r=T7N  4 

iU  READ  (1,3  )  F3  (I  +  N)  ,Z»Z1,Z2  I 

’  DO  2J  1=£73  "  | 

•20  READd,*)  Z,Z1,F2  (I)  ,Z2,Z3,Z*  | 

~3  torfitt  rrciim  * 

4  FORMAT  (SG12. 4)  • 

“I - PRinT+TTCTfr^E  l  T  0- GT3njR-TYTrE-lTTCr  "$'T'0"P . . 

KEAO*,IFG 

- rm'FG'TcQTJTSTOf5 - 

CALL  INI T 

ETto  i 

END 

snaworiNE  urn  " 

C  INITIALIZATION  ROUTINE 

~c - KEA'D"iT>r“orfA'To_TE“nrrtD - 

C  MAKE  INITIAL  GUESS 

roa  FTCfliTV'A  $/TIWuUfi  2, EE,  lP|H,N,bEL,Nl,N2  } 

COMMON/M  AT /F2  (13  2  4) ,A(2) ,AA(2> ,8(1024,2) ,D(2) ,F3(4Q96>  1 

P RfN f *73 0  FTTyPI E  TFT'  L ftUTT I m"2  ,  E E  ... 

READ*,  LI  Ml,  LI  M2, EE 

PRTn'T*  ,3H  N=,N,2HH=,H,5HLIMi=,LIhl,!3HLlM2=,LiM2,3HEE=,EE,3HIP=7: 
p:RINT*,30H  TYPE  IN  INITIAL  GUESS  -  KGAIN,  GAMMA1 
”  REAO*,A4 

PRINT*, 19H  1ST  GUESS, AA 

Bel~=i' 

CALL  COEF  (AA,  RJ2  ) 

P'RTnT  ^TTFTTmI  n  l  = ,  r 

CALL  LOO F  (RJ2 ) 

- tw- - - - - 

SUBROUTINE  DO 

~c - cAL'jerrA  rE“j  mru  r  Tnm^rmnt - 

OIMENSI3N  R (2 ) ,0 ( 2 , 2) , W1 ( 2, 2 ) 

u'C5M^WA'R/[:rH'iTm'27TE7I'PS"i7Trr0Er,T117Fre - 

COHMON/M  AT/F2 ( 1 J  2  4 ) , A (2 ) ,  AA ( 2) , 3 { 1G24,  2) ,0(2) ,F3 ( 4  096) 

- - - 

NCIM=IPjNDIMl  =  I°  +  l 

- ciQ“f3~PiTrp - 

R  (I )  =0  • 

- D*0~ ICTT^ITI'P - 

10  Q  (I ,  J )  =J  » 

- D'0~3r'K=7Tr7tT2 - 

SN=F3(K)  ~F3  (K  +  N) 

- D0-3T-Tr-lTIp - 

00  25  J=1,IP 

~?5 — (rrr77)^rr,7rmTK7*ii^5TK7'Ji - 

30  R  (I )  =R(I )  +  3  (K » I)  *  SN 

- cwnruE - 

CALL  GMINV(IP,I3,Q,W1,MR,0) 

- IT1-^rp<Tp - 

■r*IS  FAftK  IS  BIST  QOALm  W 

63  Qiii'i  f-ASl-StXSC  so  EfflO 


S 


QO  50  1=1, IP 

- 11=! - : - 

0  Cl )  =  0 . 

- DTTS3  "J=r,TP17'IP - 

0  (I )  =  0  (I )  +W1(J)*R(II) 

tj-5 - II-II+T - - 

50  0  Cl ) =-0 ( I ) 

- Eflu - 

SUBROUTINE  L00P(RJ2) 

C  ITERATION  PROCESS 

C _ COMPUTE  A  ( 1+1 )  =A  Cl)  -*-0(1) _ 

C  DETERMINE  'WHEN  A(I  +  1)  IS  ACCEPTABLE 

DIMENSION  000(2) 

COMMON/V AR/LIH1,  LIM2 ,EE , IP,H , N, DEL, Nl, N2 

COMMON/M  6T/F2  ( 1 J  24 )  ,  A  (2)  ,  A  A  ( 2  >. ,  B  ( 10  24, 2)  ,0(2) _ 

NCT  =MCT=  0 . 

1  CALL  DO _ ; _ 

2  DO  100  1=1, IP 

_ 0 ja  )  =OtI)»OEL  • _ 

1U0  A  (I)'=AA(I)+Q(I) 

IF(A(1)«  L7.6. » OR . A ( 2  ? . GT  «  0 , ) GO  TO  15 _ 

CALL  COEF  (A,'RJ1) 

6  1F(RJ1.LT. 2J2)G0  TO  30 

15  DEL=0EL/2 

_ N  CT =NCT+  1 _ _ _ 

IFtNCT.LE.LIMDGO  TO  2 

PRINT*, 7HERR0R  L,3H  A  =  ,A,5H  NCT  =  ,NCT,5H  MCT  =  , MCT , RJ1 ,RJ2 
RETURN 

_3_U _ DO  35  1=1,  IP _ 

DOO  (I )  =A  bsTo  fl )  /  A  A  ( I ) ) 

35  IF ( DOD (I ) .GT .EE)  GO  TO  20 

GO  TO  .4) 

20  MCT  =MCT+ 1 

I F ( MCT  * L  E  *  LI M 2 )*j  0  TO  9  ” 

PRINT*, 7HERROR  2,3H  A=,A,5H  NCT=,NCT,5H  MCT=, MCT , RJ1 ,RJ2 
RETURN 

j _ r  j2=rji _ ; _ 

NCT  =  d 

_ GEL  =  1 . _ 

DO~  2  5  1=1, IP 

25  A A ( I ) = A ( I ) _ 

GO  TO’  1 

40  PRINT A 5,  (A(I) ,1=1, IP) ,NCT ,MCT 

PRINT*  »3H_AA=,  A  A 
PRINT*, 3H  0= , 6 

PRINT*, 5HJMIN=TRJ1 
45  FCRHAT(iX,3G12.h,2X,2I5) 

END 

SUBROUTINE  COEF(W,RJ) 

C  COMPUTE"MEAN  TRACKING  ERROR  FOR  GIVEN  'PARAMETERS  - 

C  EXPRESS  PARTIAL  DERIVATIVES  OF  MEAN  TRACKING  ERROR  WRT 

JO  kgai'n  and  gamjai  in  "matrix  8 


64 


mis  PJVJ2  15  -v>  DOO 


DIMENSION  W(2)  ,IWK(12) 

COM  MON /VAR/ L*I  Hi»  Ll'H  2~,'EE  ,  I P  ,H  ,  N »  DEL  ,  Nl ,  N2 

_ COMMON/M  AT/F2  (1J  24  )  ♦  A  (2)  ,  AA(2 )  ,B  <1024 ,2 )  , D  (2)  ,F3  <4096) 

U  =  W  <  1 )  "  ' 

V=W{2) 

X^iJTv 
REWIND  i 
DO  10~T=T71 
WRITE(l)  F3  (I  +  N) 

T  s< 1-1 ) *  H 
B  (I )  =  F2(  I ) 

S  f=  S  2=  1-3  0T 

lF<un.GE,2‘Ji).  )31  =  0. 

"WT .L E . - 2il 0 . )  S2~=l)~ 

IFtSi.NS .C.)31=SXP(-U*T) 

IFTS2.NE.  0«)32=EX°  ( V4T) 

10  F3<I)=(tX*T+l»)*Sl«S2)/<X*X) 

Cal l  v coTfv d(F3,3,‘n,n, fwto  T 

OO  13  1=1  ,N 

- FT=T3TTT*T< - 


C 

TT 


15 


20 


60 


7u 
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H=T IMF  STEP 
TNP'GTTrAjTC 
WKITE(1)FF 
B  (I )  =F2( I ) 


04 


T=(I-1)*H 

ST=S'2=TnT7 


IF<U*T,GE.2J0.)31=0. 

TfJV*  T7lTT::IFS  . )  S  2~=TT 
IFIS1.NE.  G.)Si=EXP(-U*T) 
TTrsTtNE  .  07)  S 2 -E X P~IV**f] I 
F3<I)=((X*T-i»)*S2+Si)/<X*X) 
'CAUITTC  JiTVTtr  37T7^Fr, W) - 

DO  20  1=1, M 

TT'smrrfH - 

WRITE f  1)  FF 

■BTIWrm - 


T  =  ( I  - 1 )  ♦  H 

STi-S^^TTiD"; - 

IF(U*T «GE«290.)31=0, 
THT^TTurr^TOTr^ir; — 

IFfSl.NE. Q,)Si=EXP  <-U*T) 

TFTsTT  i;TzT(r;TST=  il'PTv  ♦  TT"  ' 

F3( I )  =  (3  2-31)  /X 

rarrT^wjTF  3iT,’NTf^7im 

REWINO  1 

'STDsir; - 

DO  60  1=1, N 

■REATDTDT  TTT+nI - 

CONTINUE 

W73”Wi;MZ - 

F3(  I) =F3 ( I)  *H 

i<j^Rj+ rn  (t+n  )  -r  3  nrn  **z 

OO  65  1=1,2 
OO  65  J=  1 , N 
R  EA  O  ( 1 )  3  ( J ,  I ) 

TO - - - 


is  B£ST  QUALIFY 
mm  i&c 


Comments:  Execution  of  Program  Fit  1. 


1.  Compile  Program 

2.  Input:  Tape  1  contains 

a.  Azimuth  and  elevation  trajectory  information 

b.  Empirical  standard  deviation  of  tracking  error 

3.  Load  Library  Routines: 

a.  VCONVO  -  convolution  integral  computation  from  International 
Mathematical  &  Statistical  Library  (IMSL). 

b.  GMINV  -  matrix  inverse  routine  from  D.  L.  Kleinman  Library 

[12]. 

4 .  Execute 

5.  Output:  Parameter  values  for  [a  ,  a  ,  a  ] . 

12  3 


67 


_ _.B.R.05.RA  M_£ITi(  INPUT., OUTPUT  ,JA  P-El) _ 

C  INPUT  MEAN  TRACKING  ERROR.  PARAMETERS 

_ P.RI N.TJU3U H . JJf.RE_IALG.AMMAi_.iKG AIN _ 

READ*, Alt A2 

_ REHXNQ_i _ 

CALL  CONVO( Alt A2) 

1 _ P.RI.NI*.t5.CHI.Y.P.L_l_Tja_a0._.QR_r.y.R.E_iJ-T.0_SJi3e. _ 

READ* , IFG 

_ LE.lIE.G.tii.Q^Q.  J.SLT.QE _ 

CALL  I  NIT 

_ 3,P„J.P__1 _ 

END 

_ S..U.3R0  ILT.INE— I.N  LI _ 

C  INITIALIZATION  ROUTINE 

C _ REAIl^JLJDi\T.A..JLQ_B£_EJ.TJJID _ _ _ _ 

C  MAKE  INITIAL  GUESS 

_ C.QHMON./J/.AR/.LIMl.,.LI.M24  EE.tH4jN.-t-D  E-Lt.Nl.tN2 _ 

"  COMMON/MAT/A<3) , AA (3 ) j B (1024, 4 > , D ( 3) ,F2<1024) 

_ LR=3 _ 

PRINT*, 3CH  TYPE  IN  LlMl,LIM2,EE 

_ &  EADl.t.LlMl ,  LIH2,.E£ _ _ 

PRINT*, 3H  N=,N,5HLIHl=,LIMi,5HLIH2=,LIM2,3HEE=tEE,3HIP=, IP 

_ -E.RINT *..».S.QH_ XYE£._I N I II AJ GUESS  =jrAL P H A1 , A LPHA 2  ,ALPH A3 _ _ 

READ*, AA 

_ PRINT.*  ,  IE  H...1.ST_.£.U£SS.,AA _ 

D£L=1 

_ v-CALI CD  EEL  A  A.,  RJ2J _ „ _ 

PRINT*, 1GH  JMIN(i) =  ,RJ2 

_ C  A  LL  -LOOP-  (RJ  2.,  IPJ _ 

END 

_ SUQRO.U.IINE_.DD.tIP.I _ 

C  CALCULATE  D  MATRIX  FROM  Q  ANO  R 

_ aiHENS.lQN..R.i.3.)„,Jl.t.3.,.3.).J.H1134.31 _ _ 

COMMON/VAR/ LIM1, LIM2 ,EE ,H,N , DEL, Nl, N2 

_ D.OKMQ N /MAT /A  (XL,  AA13J..,J3.(lfl 2  44-4  }.,£.( 3J_tE.2Lia2.4J _ 

COMMON/MAINl/NDIM,NDIMi 

_ NMM.sXPLNOI  M1=.I  P.+1 _ 

DO  ID  1=1, IP 


10 

DO  10  J=1,IP 

0  (T..1)  =f>- 

DO  35  K=N1,N2 

.  SN=B(Kt41-F2MO  . 

DO  31)  1  =  1, IP 

.D.O  25 . J=1,IP . . . . 

25 

Q(I,J) =0(I,J)+8(K,I)*B<K,J) 

30 

R(I)  =.R  (I)  +B  (  K«  I)  *SN  . . 

35 

CONTINUE 

- CALL_G.MINV.t  LP_,  IR,_Q,  Wl,  MR.O ) . 

IP1=IP*IP 

_ D.Q_5Ji  1=1.  IP  _ 

11=1 

65 


jSOil  eopy  K)  I/liQ 


0(11=0.  . . . . . 

DO  45  J=I,IPi,IP 

D(I»=0<I)+Wi(JJ»R(ID  _ 

45 

JJL_ 

11=114-1 

D  (I)  =-D  ( I )  . . . _ _ 

END 

SUBROUTINE  LOOP(RJ2.  IP)  .  ..  _ .  . 

C 

„C_ _ 

C 

ITERATION  PROCESS 

COMPUTE  A (1*1) =  A  CI)  +  D(I)  . . 

DETERMINE  WHEN  A(I+1)  IS  ACCEPTABLE 

DIMENSION  ODD (3) 

COMMON /VAR/ LI Mi»LIM2»EE»H,N,DEL»Ni»N2 

C  OUMPN  /  M  A 1/  A 

NCT=MCT=0  « 

1 

CALL  DD(IP) 

2 

DO  IDO  I=i,IP 

D  (I ) =D ( I)  ♦CEL  .  . . _ . 

Too" 

A(I)=AA(IJ+D(I) 

JF(Atl>*LTf3..0R.Af2)«LT.O.  .  OR.A(.3)_..-LT_*Q».)GO  TO  7 

6 

CALL  COEF(A,RJi) 

IFtRUJ,  ,  LT-R.12)GO  TO  30  . 

7 

DEL=0EL/2 

IF(NCT.LE.LIMi)GO  TO  2 
RETURN 


- 

15-.-. 

25 

DCrD~m=A8S(0<I)/AA<I>> 

I.F(pOOm  ,pT.EE)GO  TO  2D  .  .  . 

GO  TO  40 

MCT=MCm  . 

ifimct.le.limdgo  to  9 

PRI NT*  «  7H ERROR  2.3H  A=.A,5H  NCT=«  NOT  « 5H  MCT=,MCI,RJi,RJ2 

J3 _ 

RETURN 

?J2=RJi  .  .  . . . . 

NCT=t 

DFL=1 .  . .  . 

DO  25  1*1, IP 

25 .......  kkllY*AiU _  __ 

GO  TO  1 

PRINT*. 4H  A*.A.4HNGT=.NCT.4HMCT=,MCT  .  .  . . . .  ... 

PRINT*, 3HAA=,AA 

PRINT*. 3H  0=»D _ 

PRINT*, 5HJHIN=,RJ1 

END _ 

_C 

SUBROUTINE  COEF(W,SUM) 

COMPUTE  TRACKING  ERROR  VARIANCE  EQUATION 

c 

EVALUATE  COST  FUNCTION 

DIMENSION  W ( 3) 

COMMON/ VAR/LIM1,LIM2, EE, HiN, DEL, N1,N2 

COMMON./ MAT/ A  ( 3  ,  AA  <  3 ) ,  8  <  1 02.4 ,4. 0 13 >..» f2 $.10 Z 4 >.  . . . . . . . . . . 

feg  PAJ-2  li  Bjfc.il 

T  'Wil  COPY  MJHih-u  ?y  dd-j 


69 


SUM=0. 

_ CQ_JL0_J=fU.»J12. _ 

B(I,4) =B(I,1)*W(1)+B(I,2)*W(2) +B(I,3)*WC3) 

_ 5JJM=SU^J:tejJLiAL-E2.aJ.)..t*2 _ , _ „ _ 

10  CONTINUE 

_ END _ ; _ _ _ 

SUBROUTINE  CONVO(Al,A2) 

5. _ ,XOMPUI£_PA!LUA£_pj^XV4II.yj;S.JlF„VAR^^ 

C  DERIVATIVES  ARE  CONSTANT  SINCE  GAMMAi  AND  KG AIN  ARE  KNOWN 

_C _ LNP.UJ _ ^TRAJE.C.T.QR.lLJi^.F.ORN  A,TJ[.O.N_AN.0„_E.fl£I.RI.GAi._0.AIA. _ 


DIMENSION  IWK<12),F3<2048) 

COMMON  /  V  A  R/LIM;U.L.I  M2  tl  ZjJJ.iR  ».0..E.U.N.1j  N.2. 


COMMON/MAT/A  (3)  ,AA(3)  ,B  <102  4,4)  »D (3)  ,F2(  1(524) 

DATA  N,H,N1,N2/1G24,. 04,100,1024/ 

DO  15  1=1  ,N 

R£A0(lt*)F3(I).F2(I) »Z«Z1 

T  =  <1-1 ) *H 

D  < I ) =  < <A1*EXP(A1*T)+A2*EXP<-A2*T) ) / < A1+A2 ) ) **2 

15 

CONTINUE 

CALL  VCONVO ( 8 , F3 , N, N,I WK)  . 

DO  35  1=1, N 
.  F3<I)=F2<I) 

35 

REAO(l ,*) F2<I)  ,Z 

REWIND  1 

0d~2Q  3«1»N 
BB=S<I)*H 
WRITE(l) BB 


_ i=JI-l.)_*H _ _ _ 

20  B (I) =(  < A1*EXP( Ai*T) ♦A2*EXP(-A2*T) )/< A1+A2) ) **2 

_ C-ALL_V.CQN.VO  t.B,.F_3.,il,N,.IRK) _ 

DO  25  1=1  »N 

_ T.=  CI.-U1.H _ _ _ 

BB=B(I)*H 

_ WRlTJE.Ul.aB _ , _ 

F3 (I) =1» 

25 _ 3_(l.).=..(..ljAl.f£Xe.Uim±A2±£<R(-A2Li.T.liy_ULtA2JJ.5Lt2. 

CALL  VCONVO <B,F3,N,N,IKKi 

_ REHIMLJ. _ 

DO  30  1=1, N 

_ £2_U1=£2.(  I  )_♦  F.R  (X) _ 

30  B(I)=B(I)*H 

_ D.0_4Q— U=2.»l _ 

DO  40  1=1, N 

_4IL_ _ iEAOJLLLfiJUX  J ) - 

ENO 


.  ,u  ft*v  t  ^ 1)110 - 
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Comments:  Execution  of  Program  OBS. 

1.  Compile  Program 

2.  Input:  Tape  2  contains 

a.  Azimuth  and  elevation  trajectory  information 

3.  Load  Library  Routine: 

a.  DSCRT  -  computation  of  transition  matrix  from  D.  L.  Kleinman 
Library  [12] . 

4 .  Execute 

5.  Output:  Time,  mean  tracking  error,  standard  deviation  for  azimuth 
and  elevation. 
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_ PROGRAM,  OBS(  INPUT, OUTPUT. TAPE 2) _ 

DIMENSION  At3,3)»Z(3),X(3»3),Ri(3)»Wl(3,3)»W2{3»3),P(7,2) 

_ Q.O.MMON  /MAIN  1/N,N? _ 

.DATA  P/-2 .87.,  -1. »  2. 9 4 , C  *»  »!)496».0ti24»  *  10  3 , —  3.  i'l,“i», 3.02,0.  ,  •  0  C  3 

1  *00547.  .,25 9/jDEL  tN,  TENOA, 5 4.3.45./ _ 

C  REDUCED  OROER  MOOEL  SIMULATION  FOR  S60  SYSTEM 

C _ PJ.1  iLl2.Lf.GAM.MAi> _ PA  2L».1£1.eGA.M  M  A  2 _ £_L3_.XCJ_=_k:GAIN_ 

C  P (5.IC) =ALPHAi  P(6,IC)=ALPHA2  P(7»IC) =ALPHA3 

fi _ HQOEL-JJSE.S,J _mT.LSJmJ3.EL=.,  04.  IIMf-SlfP— _ 

N1=N*N  $  N2  =N+1 

_ L 1=1  S  L  2=2. _ 

PRINTS, 5CH  TYPE  1  FOR  AZ  2  FOR  EL  3  FOR  BOTH 

_ 1EA0*j..IFG _ 

IFUFG..EQ.1)  L2=l 

_ miE.G_«Jia .  2AU=1 _ 

DO  500  IC  =  L1* L2 

_ 8£WIUQ„2 _ 

z  <l)=-i, 

_ i.(iL=LL3l=P.l3.».Xi2J. _ 

DO  10  J=1»N1 

10 _ LUAjrMJl-A.* _ 

A (1) =P (If IC) 

_ M7_L=i, _ 

A(2)="P<l»ICi*P(3tIC) 

_ A-L9.>.=A.UJ..=-PJ.3aJ:.C)_ _ 

C  COMPUTE  TRANSITION  MATRIX  -  Wi 

_ CAlL_0-S.CBU.lli  AjJ3JELiill.,_w  2x10  J _ 

C  CALCULATE  CONSTANT  MATRICES 

_ ^ _ D.  0 5jQ I.=JL»H _ , _ 

II=i 

mi  l=jl._ _ 

DO  45  J=I * Ni » N 

jRlJU:L=RUJi.±HfiJ.LlZ.aiJ _ 


45  II=II«-1 

50 _ CONUiiUf _ 

DO  60  1=1, N 

_ ZAULrll. _ 

DO  60  J=i»N 

61 _ A  (Jj  J)_=Ri  (IJ^RJJJl _ , _ 

DO  20  1=1, N 


N 

I2=Ii+N 

__2!l _ Ri < Il=W2iILi±W2.(X21 _ 

T=0. 

,_C _ I NPUT  TRAJECTORY  INFORMATION  AND  CALCULATE  NOISE  COtf ARIANC 

1  READ (2 , 3) Cl, C2,P (4,i),C3,C4,P (4,2) 

_ 3 _ F  ORMAJJ  6G  i2 , 4 )  _ 

P4=Z(2)-Z(3) +P ( 3 , IC)  *Z<1) 

P5=(P4-PP4) /DEL 


V=(P(5fICM-P(6,IC)4P4*P4+-P(7 ,  IC)  *P5*P5)  /  DEL 
PP4=P4 


DO  25  1  =  1, N 
11=1 _ 


H2(I)=3. 

DO  15  J=I.N1.N 

15 

W2(I)  =  W2(I)  ♦WKJ)^Z(il) 

11=11+1 

25 

C 

CONTINUE 

UPDATE  MEAN  ERROR  EQUATION 

DO  35  I  =  i * N 

Z(I)=W2(I)+R1(I)*P(4,IC) 


35 

CONTINUE 

C 

UPDATE  ERROR  COVARIANCE  EQUATION 

CALL  MULT(Wl,X,N,Ni,W2) 

OLQ-AiL. I_=X»iLJ _ 

40  X  ( I)  =  A ( I)  *  V+  W2  ( I ) 

_ sn=s.QR.T-{xa,  uj _ 

C  OUTPUT  MEAN  ERROR  TRACKING  ERROR-Z(i),  AND  STANDARD  DEVIATION  C- 


.C _ LRACXlN-Gu..ERRQR.r.SQlLAR£— RQO  T  J3E-.X  ti4-U 

T  =  T  +DEL 

_ LK=JLI±j.DIl1LZD£1 _ 


IF( MOD (LK»25}«£Q«0)PRINT75,T,Z(i) ,SD 

75 

FORMAT (5X.3G12.4)  ..  ..  _  . 

IFtT.GE.TENDTGO  TO  500 
_GO  TDi _ 

590 

CONTINUE 

END  _  . 

SUBROUTINE  MULT ( E , F, L, LI ,H) 

_ DIMENSION  E (LI)  .  F  ( LI G  (Q )  .  H  ( L  1)  .  _ 

C 

MATRIX  MANIPULATION  H=EFE« 

— HCLJL9 _ laljJ _ 

II=i 


.DH_1B._K=JU1. 
TEMP=0  . 


5 _ 

TEMP=T  EMP+E ( J) ♦FT  II) 
II=II>1 

in  _ 

KK=(K-i)*L+I 

G(_KK)  =  T£MP 

OO  20  1=1, L 

DO  2R  K=I.L 

TEMP=0  • 

.  IT  =  K . .  . . 

DO  15  J=I, L1,L 

TEHP=T  EMP+G(J)+E(IT)  _ 

15 

II=II+L 

KK=(K-i)*L+I 

20 

H(KK)  =  TEMP 

_ L2=L?J _ 

DO  30  1=1, L2 

L  3s T  +  1 

DO  30  J=L3,L 

_ Kl  =  ( 1-.1J  +  L+.1  ... _ _  _  _  ... 

_ 33 _ 

K2=(J-i)+L+I 
.  4fKl)  =  HtK?) 

END 
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